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^ ' Abstract 

Q ■ The large-distance behaviour of the adiabatic hybrid potentials is studied in the 

\Q ', framework of the QCD string model. The calculated spectra are shown to be the result 

^ ; of interplay between potential-type longitudinal and string-type transverse vibrations. 

> ; 

£T) [ General arguments from QCD and lattice data tell that the theory, even quenched in 

quarks, possesses nontrivial spectrum, so that effective degrees of freedom for constituent 
glue should be introduced to describe QCD in the nonperturbative region. As far as we 
know the possibility for mesons with gluonic lump to exist was first considered in JIJ in 1976. 
Modern wisdom tells that the area law asymptotics for the Wilson loop implies a kind of 
string to be developed between quark and antiquark at large distances, and it is natural 
Q_i! to identify the qq system connected by the string in its ground state with conventional qq 
meson, while the string vibrations are responsible for gluonic (hybrid) excitations. This 
picture, though physically appealing, does not follow directly from the QCD, and one relies 
upon models to describe these excitations. There are two main ideas on how to construct 
such models. One is to consider point-like gluons confined by some potential- type force [0, [| , 

^ \ and another is to introduce string phonons Q. 

In principle the best way to discriminate between these two possibilities is to compare pre- 
dictions with experimental data on hybrid mesons. Indeed, there is a lot of indications that 
hybrid mesons are already found, but the conclusive evidences have never been presented, 
nor have alternative explanations been completely excluded ||. 

On the other hand, lattice calculations are now accurate enough to provide reliable data 
on the properties of soft glue and to check the model predictions. In this regard recent 
measurements || of adiabatic hybrid potentials are of particular interest. These simulations 
measure the spectrum of glue in the presence of static quark and antiquark separated by 
some distance R. Not only these potentials enter heavy hybrid mass estimations in the Born- 
Oppenheimer approximation. The large R limit is important per se, as the formation of 
confining string is expected at large distances, and direct measurements of string fluctuations 
become available. It is our purpose to investigate the large-distance behaviour of adiabatic 
potentials in order to establish what kind of the effective string degrees of freedom are excited 
at large distances. 
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We perform these studies in the framework of the QCD string model. This model deals 
with quarks and point-like gluons propagating in the confining QCD vacuum, and is based 
on Vacuum Background Correlators method J7J. The QCD string model was successfully 
applied to conventional mesons [f|, hybrids || 0, [7], glueballs [|llj and gluelump (gluon 
bound to the static adjoint source) fT2fl . 

The QCD string model for gluons is derived from the perturbation theory in the non- 
perturbative background, developed in [O. This formalism allows to introduce constituent 
(valence) gluons as perturbations at the confining background. The latter is given by the 
set of gauge-invariant field strength correlators responsible for the area law. The main fea- 
ture of this approach is that, in contrast to the above-mentioned models, here one is able 
to distingiush clearly between confining gluonic field configurations and confined valence 
gluons. 

The starting point is the Green function for the gluon propagating in the given back- 



ground field |L3 



G^x, y) = (D 2 (B)5, l/ + 2igF fa ,(B))-\ (1) 



where covariant derivative D x a (B) is 



D c x a (B) = 6 ca d x + gf cba B{. (2) 

The term proportional to F /ll> (B) is responsible for the gluon spin interaction; in these 
first studies we neglect it, as it can be treated as perturbation |TTl [I2||. The next step is to 



use Feynman-Schwinger representation for the quark- ant iquark-gluon Green function JTO 
which is reduced in the case of static quark and antiquark to the form 

G(x g ,y g ) = J o dsj Dz g exp(-K g )(W) B , (3) 

where K g — | J S z 2 (r)dr, and all the dependence on the vacuum gluonic field B^ is contained 
in the Wilson loop 

w = MK^ q h^f g (y 9 ,x g ). (4) 

Here <& 9 and <J>q are parallel transporters 

® q = Pexpig B^Zqjdzq^ ^g = Pexpig B^z^dz^ (5) 



with integration in (5) along the classical trajectories z qfJ , = (r, ^) and z gil = (r, —5=) of 
static quark and antiquark, P means path ordering, and 

$r 6 s (%>* 9 ) = (Pexpig f B^z g )dz g ,) a \ (6) 

a, b are adjoint colour indices, A a are Gell-Mann matrices and the contour T g runs over the 
gluon trajectory z g . 

The main assumption of the QCD string model is the minimal area law for the Wilson 



loop average, which yields for the configuration (4) the form (TC 



N 2 - 1 

(W)b = — ~ — exp(-a(St + S 2 )), (7) 
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where S% and S 2 are the minimal areas inside the contours formed by quark and gluon and 
antiquark and gluon trajectories correspondingly, and a is the string tension. 

With the form (7) for (W) b the action of the system can be immediately read out of the 
representation (3): 

T ( '2 1 

A = l dT \~2 + ^Y~ a L d /WK™0 2 - 1^?- 



-a 



£ dfa^(w 2 wtf-w 2 2 w' 2 2 } , (8) 



where the minimal surface S\ and S 2 are parametrized by the coordinates w^r, fa), i 
1,2, w ilt = ^, w> - dw ^ 



In what follows the straight-line ansatz is chosen for the minimal surface: 

wu, = t, w 1)2 = ±(l-/3)- + /3r. (9) 



The quantity fi = /i(r) in the expression (8) for the action is the so-called einbein field 
here one is forced to introduce it, as it is the only way to obtain meaningful dynamics for 
the massless particle. Moreover, we introduce another set of einbein fields, V{ = Ui(r,fa) to 
get rid of Nambu-Goto square roots in (8) ||. The resulting Lagrangian takes the form 

JO Ll><y JO I 



R 



ll,2 = f 2 -—(r 1 , 2 r)\ r 1>2 = r±^. (10) 

r l,2 1 



It is clear from Eq.(10) that the einbein field fi can be treated as the kinetic energy 
of the constituent gluon, and the einbeins Vi(r, fa) describe the energy density distribution 
along the string. These quantities are not introduced by hand, but are calculated in the 
presented formalism. Indeed, as no time derivatives of the einbeins enter the Lagrangian 
(10), it describes the constrained system, with the equations of motion 

dL - 5L - (11) 

dfi ' 6vi(fa) 

playing the role of second-class constraints. 

Now one obtains the Hamiltonian H = pr — L with the result 

a, 



2(/x+ Ji + J 2 ) 
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2A(/i + J 



vP r 2, 



^J(r 1 r 2 )(pr 1 )(pr 2 ) 

' 1' 2 



" J 2 (/X + </ 2 ) + 



(13) 



A = (// + Ji)(// + J 2 ) -J x J 2 



r\r\ 



Ji 



As we deal with the constrained system, the extra variables \i and v\p should be excluded 
by means of the conditions 

^ = °- rm = ° < 14 > 

ay. ovi{pi) 

before quantization; the extrema of the einbeins should be found from the equations (14) 
and substituted into the Hamiltonian. Such procedure is hardly possible analytically with 
the complicated structure (12), (13) even at the classical level, and after quantization these 
extremal values of einbeins would become nonlinear operator functions of coordinates and 
momenta with inevitable ordering problems arising. In what follows we use the approxima- 
tion which treats y and Ui as c-number variational parameters. We find the eigenvalues of 
the Hamiltonian (12) as functions of y and Ui and minimize them with respect to einbeins 
to obtain the physical spectrum. Such einbein method works surprisingly well in the QCD 
string model calculations, with the accuracy of about 5-10% for the ground state |L5| . 

Even with this simplifying assumption the problem remains complicated due to the pres- 
ence of the terms Ji )2 responsible for the string inertia. Suppose for a moment that one can 
neglect these terms in the kinetic energy (13). Then the Hamiltonian takes the form [|TU|, U\ 



2^,2 



or 



u 2v 2 



+ 



2y 2 Jo 2v x 

which allows to eliminate einbeins and to arrive at the potential model Hamiltonian 



(15) 



H 



/p 2 + err i + or 2- 



(16) 



Let us now estimate whether the neglect of string inertia is justified. To this end we find 
the spectrum of the Hamiltonian (15), (16) using the einbein method described above. It is 
given by the set of equations 

, s Mn + I) 2 c> 2 
E n (R)=^ n (R)+ K 2) 



16o> + ^) 4 = A4(i2)(4(n + 3 -) 2 + R^ 2 n (R)) 



(17) 



with Vi independent of 



, s 2(n + |)V 2 



MR) ' 

is the projection of orbital momentum onto z axis, z || R. 
Note that while the angular momentum is not conserved in the exact Hamiltonian (16), it is a 



where n = n z + n + A, A 



LR 



R 
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good quantum number in the approximate einbein method: we have compared the spectrum 
of exact and einbein-field Hamiltonian and have found that angular momentum is conserved 
in the potential problem (16) within better than 5% accuracy. The same phenomenon 
is observed in the constituent gluon model |§, and is the consequence of linear potential 
confinement. 

Consider first the small R,R^ l/y/a, limit of the system (17): 

gn (ii) = 2^(n + |)^ + 23/2 ^ )i/2 , (19) 

o 3/2 d2 

/!„(#) = 2 1 /V/ 2 (n + -) 1 /2 



2 J (n + |) 1 /2 2 5/2 



, D . (n+D^V/ 2 3a 3 / 2 R 2 

v\,2n(R) = T^n 1- 



2V2 27/2(^ + 1)1/2- 

The last line in (19) yields Jip/ A 4 ~ g- ^ ne situation here is similar to the one in the light 
quark, glueball and gluelump QCD string calculations: the correction due to string inertia 
is sizeable but not large, and can be taken into account as perturbation JTI, T^. Note that 
it is the regime of small R which is relevant to the heavy hybrid mass estimations fl6|| : the 
average distance between heavy quark and antiquark is small, (R 2 ) <C 1/cr, so that QQ pair 
resides in the oscillator adiabatic potential which, in the einbein method, is given by Eq. 
(19). 

The situation changes drastically for the case of large R, R^> \j \fo. Now gluon enjoys 
small oscillation motion, and one has 

MR) = °R+iry l3( ^^, (20) 

4a'>. + j)^ oR 

^n{R) - -^pi , V\,2n\R) - — , 

displaying (J^) 1 / 3 subleading behaviour typical for linear potential confinement at large dis- 
tances j7|. Nevertheless, in this case Ji )2 = \&R 3> Hn, so that the potential regime is 
unadequate at large R. 

To get more insight into what happens at the intermediate and large distances we consider 
the quasiclassical limit of large A, where only rotations around z axis are taken into account: 

H _ A 2 , 
2p 2 (J 1 + J 2 ) 

r^T+r^T- < 2i) 

As no momenta p z and ^j- enter the Hamiltonian, the system stabilizes itself at the points 
z and po given by the conditions 
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Combining Eq.(22) with the second condition of Eq.(14) one arrives at the following expres- 
sions: 

,20 = 0, /0o = — — , v L (/3) = u 2 (J3) = v(J3), (23) 
lay J a 



where 

J= f 1 d(3(3 2 u{(3), a= f 1 —j-^r (24) 



o 



and the function v{(3) is given by 



VA a 2 R 2 Act Act /a" , . 

Substituting the form (25) into eqs. (24) one finds the expression for the energy 



,„ ,/, ^{arcsinv / 5 + x /5(l-5)} 1 / 4 

E = 2CT 1/2 A 1/2 arcsin\/B- V (26) 

{arcsinv^B - J B(l - 5)} 3 / 4 



2Afi 3/ Vl - 5 

ct# 2 



{arcsin\/£? + - £)} 1/2 {arcsiir/B - /B(l - £)} 



13/2 



4 

with the large R limit of (26) given by 

E(R) = aR+2y/3^. (27) 
R 

Here we have 1/R subleading behaviour typical for naive Nambu-Goto string models. For 
example, the flux-tube model j| predicts 

E(R) = a R + ^- (28) 
it 

in the small-oscillation approximation. The energy curve (26) is shown at Fig. 1 together 
with the flux-tube (28) and potential-regime curve (17) for n z = n p = and A = 1,2,3. 
The large R limit of the quasiclassical regime (26) is very close to the flux tube one and 
deviates substantially from the potential regime, while at small R unphysical divergent 1/R 
behaviour is absent. 

The case of large R can be treated directly in the full Hamiltonian (12), which in the 
small-oscillation limit takes the form 

H = ^ + ^+ Pl 



2 2/x 2(/i + Jx + J 2 ) 



+ /' dp!^ + /' *3 2 ^, (29) 
jo 2 Jo 2 
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E, GeV 




Figure 1: Adiabatic hybrid potentials in various regimes. Quasiclassical (solid line), potential 
(dashed), and flux-tube (dotted) curves for n 2 
aR. a = 0.22GeV 2 . 



n p = and A = 1,2,3. The lowest curve is 



E n -aR, GeV 
1 . 2 

C . 8 * 

. 4 



A=l 



E n -aR,GeV 



E n -aR,GeV 



R, fm 



A=5 



R, fm 



A=10 



R, fm 



Figure 2: Corrections to linear behaviour of potentials. QCD string (solid line), potential 
(dashed), and flux-tube (dotted) curves; n z = n p = 0; a = 0.22GeV 2 . 



displaying two different kinds of string excitations, along the z axis and in the transverse 
direction. Indeed, for large R one neglects the contribution of fi in the third term of (29) 
because the extremal values of are Then the oscillations in the longitudinal and 
transverse directions become uncoupled, and one has 

E n (R) =aR+ ^ [ Rl j 3 2) + =-|— K + A + 1). (30) 

The regime ~ (^) 1//3 is established at large R, but at the intermediate distances there are 
sizeble corrections from the string regime ~ A/i?, as it is seen from Fig. 2. 

As we have not considered the spin of the gluon, we are not in the position yet to compare 
our predictions with lattice results 0. Nevertheless, some preliminary conclusions can be 
drawn. For separations less than 2 fm the measured energies || lie much below Nambu-Goto 
curves (28). There is no universal Nambu-Goto behaviour even for R as large as 4 fm. The 
QCD string model is able to describe both these features: at small separations the potential 
confinement regime dominates, while at large distances the situation is more complicated. 
Indeed, there is the contribution of the string-type gaps (27) which are due to transverse 
vibrations of the string, but the dominant subleading behaviour is defined by potential-type 
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longitudinal motion. In particular, even for quasiclassically large values of A there exists the 
contribution of oscillations in the longitudinal direction (second term in (30)). 

Such peculiar behaviour displays the most pronounced difference between the given ap- 
proach and other models of constituent glue. In contrast to phonon-type models, the QCD 
string vibrations are caused by point-like valence gluon, but, in contrast to potential models, 
the confining force follows from minimal area law, giving rise, at large distances, both to 
longitudinal vibrations with potential-type ~ (^) 1//3 dominant subleading behaviour and to 
the transverse vibrations with string-type ~ A/R subleading behaviour, which could be re- 
sponsible for the observed A dependence. The full QCD string calculations with gluon spin 
involved will provide, if confirmed by the lattice data, the decisive evidence in favour of the 
QCD string model of valence glue. 

We are grateful to Yu.A.Simonov for useful discussions. The support of INTAS-RFFI 
97-0232, RFFI 00-02-17836 and 00-15-96786 grants is acknowledged. 
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